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Abstract 

A shear loaded, stringer reinforced composite 
panel is analyzed to evaluate the fidelity of 
computational fracture mechanics analyses of 
complex structures. Shear loading causes the panel 
to buckle. The resulting out-of-plane deformations 
initiate skin/stringer separation at the location of an 
embedded defect. The panel and surrounding load 
fixture were modeled with shell elements. A small 
section of the stringer foot, web and noodle as well 
as the panel skin near the delamination front were 
modeled with a local 3D solid model. Across the 
width of the stringer foot, the mixed-mode strain 
energy release rates were calculated using the 
virtual crack closure technique. A failure index was 
calculated by correlating the results with a mixed- 
mode failure criterion of the graphite/ epoxy material. 
The objective was to study the effect of the fidelity of 
the local 3D finite element model on the computed 
mixed-mode strain energy release rates and the 
failure index. 


1 Introduction 

Many composite components in aerospace 
structures are made of flat or curved panels with co- 
cured or adhesively bonded frames and stiffeners. 
Recent studies focused on the investigation of the 
debonding mechanism and included testing of 
skin/stiffener panels and failure analysis using shell 
models [1, 2]. Over the last decade, a consistent 
step-wise approach has been developed which uses 
experiments to determine the failure mechanism, 
computational stress analysis to determine the 
location of first matrix cracking and computational 
fracture mechanics to investigate the potential for 
delamination growth. Testing of thin-skin stiffened 

If Resident at Durability, Damage Tolerance and Reliability Branch, 
MS188E, NASA Langley Research Center, Hampton, VA, 23681, USA. 


panels designed for use in pressurized aircraft 
fuselages has shown that bond failure at the tip of 
the frame flange is an important and very likely 
failure mode [3]. Comparatively simple, laboratory 
size coupon type specimens consisting of a stringer 
flange bonded onto a skin have been developed to 
study the fundamental mechanisms of skin/stiffener 
debonding [4]. The failure that initiates at the tip of 
the flange in these coupon-type specimens is nearly 
identical to the failure observed in the full-scale 
panels and frame pull-off specimens [4-6]. A 
methodology based on fracture mechanics [7] has 
been used successfully to investigate the onset and 
growth of delaminations in simple characterization 
specimens and laboratory-size coupon type 
specimens [5, 6]. Future acceptance of the 

methodology by industry and certification 
authorities, however, requires the successful 
demonstration of the methodology on a structural 
level. 

For the demonstration of the methodology on 
the structural level, a stringer stiffened panel, as 
shown in Fig. 1, has been analyzed previously [8, 9]. 
The square (1016 mm x 1016 mm) panel made of 
carbon/epoxy tape is reinforced with three stringers 
made of carbon/epoxy plain weave fabric. Details 
will be discussed in the following sections. During 
manufacturing, an artificial defect of about 82 mm in 
length was placed at the termination of the center 
stiffener. The stiffened panel was bolted to a steel 
picture frame and subjected to shear loading which 
caused the panel to buckle as shown in the finite 
element model depicted in Fig. 1 [8, 9]. The 
resulting out-of-plane deformation caused 
skin/stringer separation to initiate at the location of 
the artificial defect. A small section of the stringer 
foot and the panel skin in the vicinity of the 
embedded defect were modeled with a local 3D 
solid model as shown in the enlargement in Fig. 1. 
The mixed-mode strain energy release rates were 
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calculated using the virtual crack closure technique 
[10, 11] across the width of the stringer foot. A 
failure index was calculated by correlating the 
results with the mixed-mode failure criterion of the 
graphite/epoxy material [8, 9]. 


peak 


w=0 


trough 


u=v=0 


panel made of 
carbon/epoxy tape 
modeled with shell 
elements 

gap in load frame and 
slit in panel corner 


steel load frame 
z,w modeled with shell 
and beam elements 




stiffener made of carbon/ 
epoxy fabric modeled 
with shell elements 


i=v=6.35 mm 



Fig. 1. Buckling of stringer stiffened panel subjected 
to shear loading (1016 mmx 1016 mm) 


The objective of the current research was to 
study the effect of the fidelity of the local 3D finite 
element model on the computed mixed-mode strain 
energy release rates and the failure index. In the 
original model, shell elements represented the skin 
and stiffener foot [8, 9]. In previous studies, a model 
was introduced where only the stiffener foot and 
skin were included in a local 3D model, while the 
web and hat were modeled with shell elements as 
shown in Fig. 1. The pair of 2D illustrations in Fig. 2 
are cross sections of the local 3D finite element 
model [8, 9]. For the current study, a new set of 
models was introduced, as shown in Fig. 3, where 
the stringer web was included in the local solid 
model. Another set of models included the transition 
radius between the web and foot as well as the 
detailed noodle region in the local solid model as 
shown in Fig. 4. Details of the models are discussed 
later. Thus, this study considers the effects of 
various geometrical simplifications that might be 
considered by an analyst when modeling debonding 
in complex structures. The work discussed herein is 
intended to be used as a guide to determine the trade 
off between modeling complexity and fidelity on the 


corresponding results. The current research 
complements previous studies [8, 9]. 
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Fig. 2. Local 3D model of panel skin and 
stiffener foot 



Model with refined mesh near the 
center (TFL3-20/20/20) 
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solid elements 


stiffener foot and skin 
modeled with solid elements 


plane of 
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Fig. 3. Local 3D model of panel skin, stiffener foot 
and web 
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Model with refined mesh near the 
edges (TNI L3C-8/1 0/80) 


stiffener web, 
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with solid elements y> v 



stiffener foot and skin 
modeled with solid 
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plane of 
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Fig. 4 Local 3D model of panel skin, stiffener foot 
and web including the noodle and transition radius 


2 Methodology 

2.1 Linear Elastic Fracture Mechanics 

Linear elastic fracture mechanics has proven 
useful for characterizing the onset and growth of 
delamination in composite laminates [7] . When 
using fracture mechanics, the total strain energy 
release rate, G T , is calculated along the delamination 
front. The term, G T , consists of three individual 
components. The first component, G h arises due to 
interlaminar tension. The second component, G n , 
arises due to interlaminar sliding shear (shear 
stresses parallel to the plane of delamination and 
perpendicular to the delamination front). The third 
component, Gnu arises due to interlaminar scissoring 
shear (shear stresses parallel to the plane of 
delamination and parallel to the delamination front). 
The calculated G h G n , and G IU components are then 
compared to interlaminar fracture toughness values 
in order to predict delamination onset or growth. The 
interlaminar fracture toughness values are 
determined experimentally over a range of mixed- 


mode ratios from pure mode I loading to pure mode 
II loading [12-14]. 

A quasi static mixed-mode fracture criterion is 
determined by plotting the interlaminar fracture 
toughness, G c , versus the mixed-mode ratio, G n /G T 
as shown in Fig. 5. The fracture criteria is generated 
experimentally using pure Mode I (G n /Gf= 0) 
Double Cantilever Beam (DCB) tests, pure Mode II 
(G n /Gr= 1) four point End-Notched Flexure (4ENF) 
tests, and Mixed Mode Bending (MMB) tests of 
varying ratios of Gj and G n . A fracture criterion was 
suggested by Benzeggah and Kenane [15] using a 
simple mathematical relationship between G c and 
Gjj/Gt 


G,-G k + (G„,-G„y 


(G V 

—A 

\G T ) 


( 1 ) 


In this expression G Ic and G IIc are the 
experimentally-determined fracture toughness data 
for mode I and II. The factor 77 is determined by a 
curve fit. Fracture initiation is expected when, for a 
given mixed mode ratio G n IG T , the calculated total 
energy release rate, G T , exceeds the interlaminar 
fracture toughness, G c . 



Mixed-Mode Ratio G || /G J 

Fig. 5. Mixed-mode failure criterion 

Several specimens, including the edge-cracked 
torsion specimen (ECT), have been suggested for the 
measurement of the mode III interlaminar fracture 
toughness property [16, 17]. However, mixed-mode 
I/III and II/III tests are as of yet unavailable, and 
thus a complete 3 -dimensional interaction criterion 
such as that suggested in [18] cannot be accurately 
defined. Therefore, a modified definition is 
introduced for three-dimensional analysis, which 
also yields results for the scissoring mode G m . In the 
modified definition G s denotes the sum of the in- 
plane shearing components G n +G n i [8, 9]. This is 
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necessary since a mixed-mode failure criterion, 
which accounts for all three modes, is currently not 
available. For analyses where G n f= 0, this definition 
is equal to the commonly used definition of the 
mixed mode ratio, G n /G T mentioned above. 

To determine failure along the delamination 
front, the critical energy release rate G c is calculated 
using Eq. (1) with G n = G s at each point along the 
delamination front. Subsequently, the failure index 
G t /G c is determined from the computed total energy 
release rate, G T , and the critical energy release rate 
G c with the assumption that delamination 
propagation occurs for 

G r 

— ^ 1 . ( 2 ) 
G c 

2.2 Analysis Tools 

2.2.1 Virtual Crack Closure Technique 

A variety of methods are used in the literature 
to compute the strain energy release rate based on 
results obtained from finite element analysis. For 
delaminations in laminated composite materials 
where the failure criterion is highly dependent on the 
mixed-mode ratio, the virtual crack closure 
technique (VCCT) [10, 11] has been most widely 
used for computing energy release rates. VCCT 
calculations using continuum (2D) and solid (3D) 
finite element analyses provide the mode separation 
required when using the mixed-mode fracture 
criterion. 

2.2.2 A Global/Local Shell 3D Modeling Technique 

Built-up structures are traditionally modeled 
and analyzed using plate or shell finite elements, as 
shown in Fig. 1, to keep the modeling and 
computational effort affordable. The fidelity of the 
resulting computed mixed mode strain energy 
release rate components, however, depend on many 
variables. These variables include element type, 
order of the shape functions and shear deformation 
assumptions, kinematic constraints in the 
neighborhood of the delamination front, and 
continuity of material properties and section 
stiffnesses in the vicinity of the debond when 
delaminations or debonds are modeled with plate or 
shell finite elements [19]. These problems may be 
avoided by using three-dimensional models. 
However, since many layers of brick elements 
through the thickness are often necessary to model 
the individual plies, the size of finite element models 


required for accurate analyses may become 
prohibitively large. 

For detailed modeling and analysis of the 
delaminations, the shell/3D modeling technique will 
reduce the modeling time compared to that required 
to run a fully three-dimensional finite element 
model. The technique will also reduce computational 
time because only a relatively small section of the 
mesh needs to be modeled with solid elements, 
minimizing the overall size of the model. The 
technique combines the accuracy of the full three- 
dimensional solution with the computational 
efficiency of a plate or shell finite element model. 
The technique has been demonstrated for various 
applications such as fracture toughness 
characterization specimens [20], on the coupon level 
for the skin/stringer separation specimen [21] and in 
related studies for skin/stringer separation [8, 9]. The 
enlargement in Fig. 1 illustrates the regions within 
the stringer stiffened panel that are modeled with 
shell and solid elements. 

3 Finite Element Modeling 

In the current study, a finite element analysis of 
the three-stringer panel shown in Fig. 1 was 
conducted. The load frame and the three-stringer 
panel were modeled with beam and shell elements. 
A small section of the stringer and the panel skin in 
the vicinity of the embedded defect was modeled 
with a local 3D model. For all analyses, the stiffener 
hat was modeled with shell elements. 

3.1 Global Shell Model of Stringer Stiffened 
Panel 

The global model included the steel load frame 
and attachments, the panel made of graphite/epoxy 
prepreg tape, and the stringers made of 

graphite/epoxy fabric, as shown in Fig. 1. The outer 
steel load frame and the attachment bolts were 
modeled with beam elements (ABAQUS® element 
type B21) [22]. The inner steel load frame, which 
overlaps the panel edge, was modeled with shell 
elements (ABAQUS® element type S4). The 
stiffener components, such as the foot, web and hat 
were also modeled with shell elements. 

The panel skin and the stiffener foot are 
modeled as separate entities. The S4 shell elements 
are located at the panel skin and stiffener foot 
respective mid-planes. The shell elements are 
connected by beam elements (ABAQUS® element 
type B31 modeled as steel) to provide point-to-point 
constraints between the two surfaces [23]. In the 
section containing the artificial defect, the beam 
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elements were replaced by gap elements 
(ABAQUS® element type GAPUNI). The gap 
elements allow the modeling of contact between two 
nodes. The nodes can be in contact (gap closed), 
which prevents element interpenetration or separated 
(gap open) which allows the skin/stringer separation 
[ 22 ]. 

3.2 Local 3D Model of the Stringer Foot and 
Panel Skin 

The local 3D model of the stringer foot and 
panel skin was generated using solid elements 
(ABAQUS® element type C3D8I) and consisted of 
an intact section and a delaminated section with a 
fine mesh around the delamination front. Two 
examples with different mesh refinement (BF- 
8/10/8, BF-20/20/20) are shown in Fig. 2. The 
artificial defect is located at the bondline between 
the stringer foot and the panel as shown in the 
enlargements of Fig. 2. This defect was treated as a 
delamination and modeled as a discrete discontinuity 
using two unconnected nodes with identical 
coordinates on each side of the delamination. 
Contact was modeled between the delaminated 
surfaces to avoid interpenetration during the 
analysis. 

Four elements through the thickness were used 
to model the foot of the stiffener made of 
carbon/epoxy fabric as shown in the 2D cross 
sections in Fig. 2. The skin ply made from 
carbon/epoxy tape which is adjacent to the plane of 
the delamination was modeled with one element. 
The remaining 10 plies of carbon/epoxy tape were 
modeled with three elements through the thickness 
as shown in Fig. 2. 

3.3 Local 3D Model of the Panel Skin, Stringer 
Foot and Web 

The modeling fidelity was increased by 
including the stiffener web in the local 3D model as 
shown in Fig. 3 (TFL3-8/10/8, TFL3 -20/20/20). The 
local 3D model shown consisted of two separately 
meshed sections: the stiffener web and the panel 
skin/stiffener foot. The web was modeled with eight 
solid elements (ABAQUS® element type C3D8I) 
through the thickness. The two models are 
connected to each other using multi-point constraints 
(ABAQUS® MPC option *TIE). The multi-point 
constraints are used to tie two surfaces together for 
the duration of a simulation, which makes the 
translational and rotational motion equal for a pair of 
surfaces. Nodes on the slave surface (defined by the 
user) are constrained to have the same motion as the 


point on the master surface (defined by the user) to 
which it is closest [22]. 

3.4 Local 3D Model of the Panel Skin, Stringer 
Foot and Web Including the Noodle and 
Transition Radius 

The modeling fidelity was increased further by 
including the noodle region and transition radius in 
the local 3D model as shown in Fig. 4. The 
triangular region underneath the web, where the T- 
stiffener connects to the skin, is usually referred to 
as a noodle. As above, the local 3D model shown 
consisted of two separately meshed sections which 
are connected with multi -point constraints: The 
model of the T-stiffener including the web, noodle 
region and transition radius and the model of the 
panel skin and delaminated interface is shown in the 
cross section of Fig. 4. The two solid models were 
joined using the *TIE option in ABAQUS® [22]. 

The solid model of the stringer foot and 
delamination consisted of brick elements 
(ABAQUS® element type C3D8I) and included an 
intact section and a delaminated section with a fine 
mesh around the delamination front, as discussed 
above. 

Solid models representing the stiffener foot, 
the web and the noodle region consisted of eight- 
node C3D8I brick elements, with a small number of 
triangular prism elements C3D6 used to model part 
of the noodle region as shown in the cross sections 
of Fig. 4. Since the exact transition radius was not 
known, the influence of detailed local 3D modeling 
on computed strain energy release rates was studied. 
In this case, a number of models were generated 
with different foot/web transition radii. Models were 
generated for radius values of 0.254 mm (TN3L3C- 
8/10/8, TN3L3C-20/20/20), 0.711 mm (TN1L3C- 
8/10/8, TN 1 L3 C-20/20/20), 2.54 mm (TN2L3C- 
8/10/8, TN2L3C-20/20/20) and 5.08 mm (TN4L3C- 
8/10/8, TN4L3 C-20/20/20). Models with radii 2.54 
mm and 0.711 mm were thought to correspond to 
values used in manufactured panels. The small 
radius 0.254 mm (cross section shown in Fig. 4) was 
selected to determine if the computed results were 
similar to those obtained from the model discussed 
in section 3.3 which did not included the radius and 
noodle. The larger radius 5.08 mm (cross section 
shown in Fig. 4) was chosen as an upper limit. It did 
not appear to be meaningful from a design 
standpoint to assume larger radii. 
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3.5 Finite Element Model Assembly, Load and 
Boundary Conditions 

The local 3D models were joined with the shell 
model to complete the assembled model as shown in 
the enlargement of Fig. 1. At the boundaries, the 
shell edges were used to connect the shell model 
with the local 3D solid models. The connection was 
accomplished using the shell to solid coupling 
option in ABAQUS®, which allows the connection 
between non-conforming shell and solid models. 
The coupling option uses a set of internally defined 
distributing coupling constraints to couple the 
motion of a row of nodes along the edge of the shell 
model to the motion of a set of nodes defined on a 
surface of the solid model [22]. 

For modeling the experiment, which was 
performed under displacement control, uniform 
displacements u, v were applied at one corner node 
to introduce shear as shown in Fig. 1. The in-plane 
displacements u, v were suppressed at the diagonally 
opposite corner, and the out-of-plane displacements 
w were suppressed along all four edges across the 
entire width of the inner and outer steel load frame. 

3.6 Analysis Overview 

A total of 41 different model combinations 
were analyzed. An overview of all models generated 
is given in reference [24]. The following cases are 
discussed here: 

• Models were analyzed where only the stiffener 
foot and skin were included in a local 3D model, 
while the web and hat were modeled with shell 
elements as shown in Fig. 2. 

• A set of models were created, as shown in 
Fig. 3, where the stringer web was included in 
the local solid model. 

• Another set of models included the transition 
radius between the web and foot as well as the 
detailed noodle region in the local solid model 
as shown in Fig. 4. The stiffener hat, modeled 
with shell elements, was kept unchanged. Four 
different transition radii were studied 
(r=0.254 mm, 0.71 1 mm, 2.54 mm, 5.08 mm). 

• In order to capture the local failure near the 
edges, models with a locally refined fine mesh 
were chosen as shown in the cross sections of 
Figs. 2 to 4. In order to capture the local failure 
index distribution in the vicinity of the web 
termination, models with a fine mesh in the 
center were chosen (also shown in Figs. 2 to 4). 


4 Analysis Results 

4.1 Model Deformation 

Under the applied shear loading, the analysis 
predicts the buckling deformation shown in Fig. 1. 
For the simulated delamination length (a= 81.9 mm), 
three peaks and one trough can be observed in the 
panel bays adjacent to the center stiffener as shown 
in Fig. 1. More details may be found in a previous 
study [8]. 

Details of the deformed finite element models 
are shown in the enlargements of Figs. 2 to 4 after 
the entire external displacement u=v=635 mm had 
been applied. Mode I opening was observed across 
the entire width of the stringer over the entire 
delaminated length. 

4.2 Comparison of Failure Indices Computed 
from Different Local 3D Models 

For each nodal point along the delamination 
front, the critical energy release rate, G c , was 
calculated from a mixed mode failure criterion 
(Eq. 1) for the computed mixed-mode ratio, G S /G T . 
Subsequently, the failure index G T /G C was calculated 
from Eq. 2. The failure index was calculated for the 
final load increment and plotted versus the location 5 
across the width of the stringer, b as shown in 
Figs. 1 and 2. 

The computed failure index distributions across 
the width of the stringer obtained from different 
models were plotted in Figs. 6 to 12 for comparison 
with reference results. The results from the models 
with refined edges and center region (Fig. 2) were 
used as reference results in all the figures and are 
depicted as blue circles. For all cases shown, the 
failure index peaked at the edges (5=0.0 and 5=1.0) 
with an additional peak around the center (5-0.5) 
underneath the stringer web. 



location along delamination front, s 


Fig. 6. Computed failure index - 3D web results 
versus reference solution 
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First, the results obtained from a model that 
included the panel skin, stiffener foot and the web in 
the local solid model as shown in Fig. 3 were plotted 
in Fig. 6 for comparison with the reference result. 
Qualitatively, both distributions followed the same 
trend. In two areas to the left and right of the web 
(0.0 < s< 0.4 and 0.7 < s < 1.0), the results are in 
good agreement. Locally, near the web (0.4 < s< 0.7) 
the distributions differ. The peak values computed 
for the local solid model that included the web (red 
squares) exceed the reference values (blue circles) 
by about 49%. The peak location is also offset. 

The results obtained for a transition radius 
7=0.711 mm are plotted in Fig. 7. In the area to the 
left of the web (0.0 < s< 0.4), the results are in good 
agreement. In the areas to the right of the web (0.7 < 
s < 1.0), the results are higher for the models which 
included the web and the noodle (green diamonds) 
than the results obtained with the reference solution. 
An additional small plateau was observed for 
0.47 < 5 < 0.49 (G t /G c ~ 29) which was not observed 
in the results used as reference. Also, the peak 
values computed for the local solid model that 
included the web and the noodle exceed the 
reference values (blue circles) by about 67%. As 
before the peak location is offset compared to the 
reference results. 


the reference values (blue circles) by about 51%. For 
this transition radius, the peak location is almost 
identical to the peak observed for the reference 
results 

The results obtained for a smaller transition 
radius (t= 0.254 mm as shown in Fig. 4) are plotted 
in Fig. 9. In two areas to the left and right of the web 
(0.0 < s< 0.4 and 0.7 < s < 1.0), the results are in 
good agreement. Locally, near the web (0.4 < s< 0.7) 
the distributions differ. The peak values computed 
for the local solid model that included the web and 
noodle (black triangles) exceed the reference values 
(blue circles) by about 51%. The peak location is 
also offset. 



0.4 


0.6 


location along delamination front, s 



location along delamination front, s 


Fig. 7. Computed failure index - results for transition 
radius t= 0.7 1 1 mm versus reference solution 


Fig. 8. Computed failure index - results for transition 
radius 7=2.54 mm versus reference solution 



location along delamination front, s 


For the results computed for the larger 
transition radius (t= 2.54 mm), an additional local 
maximum was observed for 0.44 < 5 < 0.46 
(G t /G c ~ 15) as shown in Fig. 8. In the areas to the 
right of the web (0.7 < s < 1.0), the results are higher 
for the models, which included the web and the 
noodle (orange triangle) than the reference solution. 
Also, the peak values computed for the local solid 
model that included the web and the noodle exceed 


Fig. 9. Computed failure index - results for transition 
radius 7=0.254 mm versus reference solution 

The results computed for a larger transition 
radius (t= 5.08 mm as shown in Fig. 4) are plotted in 
Fig. 10. In the area to the left of the web (0.0 < s< 
0.4), the results are in good agreement. For the local 
solid model that included the web and noodle (red 
squares) an additional local maximum was observed 
for 5 ~ 0.4 (G t /G c ~ 5) which was not observed in the 
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results used as reference (blue circles). In the areas 
on the right of the web (0.7 < s < 1.0), the results are 
higher for the models which included the web and 
the noodle than the results obtained with the 
reference solution. For this transition radius, the 
peak value is almost identical to the peak failure 
index computed for the reference results. 



x — i i — , — , — , — i-L-i — i — i — i — i — i — i — 

0 0.2 0.4 0.6 0.8 1 


location along delamination front, s 


Fig. 10. Computed failure index - results for radius 
r=5.08 mm versus reference solution 

For closer comparison, the computed failure 
index distributions across the width of the stringer 
are plotted in Figs. 11 and 12 for all models 
discussed. In the area to the left of the web 
(0.0 < s< 0.4), the failure index is low and all results 
are generally in good agreement. The peak values 
computed for the local solid models that included 
the web and the models that included the web and 
the noodle exceed the reference values (blue circles), 
which were obtained from models where the web 
had been modeled with shell elements. The location 
where the peak failure index was observed shifted 
from model to model. The results obtained from 
models where only the web was included in the local 
solid model (red squares) were in excellent 
agreement with the failure indices computed from 
the model with the smallest radius (r=0.254 mm, 
black triangles). In two areas to the left and right of 
the web (0.0 < s< 0.4 and 0.7 < s < 1.0), the results 
are in good agreement with the reference results. 
Locally, near the web (0.4 < s< 0.7), the 
distributions differ up to 67% as discussed above. 
For models with larger radii (r=0.711 mm, green 
diamonds; r= 2.54 mm, orange triangles and 

r=5.08 mm, red squares), an additional local maxima 
or plateau was observed for 0.4 < 5 < 0.5 as shown 
in Fig. 12. With increasing radius, the peak value 
decreases and shifts to the right of the center 
0.5 <s < 0.6. Also, in the area to the right of the web 


(0.7 < 5 < 1.0), the results are lower for the models 
with a larger transition radius. 



location along delamination front, s 


Fig. 11. Comparison of results for 3D web, radius 
r=0.254mm and reference solution 



o 

BF-8/1 0/8 + BF-20/20/20 
(Figure 2) 


o 

TNI L3C-8/10/8 + 

" 


TNI L3C-20/20/20 
(Figure 4) 


A 

TN2L3C-8/1 0/8 + 
TN2L3C-20/20/20 

_ 


(Figure 4) 


□ 

TN4L3C-8/1 0/8 + 
TN4L3C-20/20/20 
(Figure 4) 


80 - 


0 Tln a aniBnaanHB 



r=0.71 1 mm 
r=2.54 mm 
r=5.08 mm 



0.2 0.4 0.6 0.8 

location along delamination front, s 


Fig. 12. Comparison of results for different radii and 
reference solution 


5 Summary and Concluding Remarks 

The skin/stringer separation of a 
graphite/epoxy composite panel reinforced with 
three stringers and subjected to shear loading was 
studied using computational fracture mechanics 
analysis. The shear loading causes the panel to 
buckle and the resulting out-of-plane deformation 
initiates skin/stringer separation at the location of an 
embedded defect. The panel and surrounding load 
fixture were modeled with shell elements. 

A small section of the stringer foot, web and 
noodle as well as the panel skin in the vicinity of the 
delamination front were modeled with local 3D solid 
models. The mixed-mode strain energy release rates 
were calculated along a straight delamination front 
across the width of the stringer foot using the virtual 
crack closure technique. A failure index was 
calculated by correlating the results with a mixed- 
mode failure criterion of the graphite/epoxy material. 
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Computed failure indices were compared to 
corresponding results where the entire web was 
modeled with shell elements and only a small 
section of the stringer foot and panel was modeled 
locally with solid elements 

The results showed the following 

• including the stiffener web in the local 3D 
model increased the computed failure index by 
about 49% compared to the reference results 
where the web was modeled with shells 

• including the web, the noodle and transition 
radius in the local 3D solid model increased the 
magnitude of the failure index up to 67% and 
changed the local distribution across the width 

• the magnitude of the failure index decreased 
with increasing transition radius. 

Differences in the failure index distributions 
were likely caused by a different local deformation 
behavior due to different local stiffnesses of the 
models studied. The local differences in stiffness 
were mainly caused by the local modeling of the 
noodle and the transition radius. Based on the 
increase in computed failure index, it is suggested to 
use a high fidelity model including the noodle and 
transition radius whenever accurate analysis results 
are required. The results of this study are intended to 
be used as a guide for conducting finite element 
analyses of structures such as stiffened panels. In 
particular, this guidance is aimed towards analyses 
that attempt to simulate delamination growth and 
debonding. 
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